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Abstract. We present results from the first generation of multi-dimensional general 
relativistic neutrino hydrodynamics simulations of core-collapse supernovae. A com- 
parison with models computed using either the purely Newtonian approximation or the 
"effective gravitational potential" approach reveals appreciable quantitative differences 
in the heating conditions and the gravitational wave spectra. Our results underscore the 
important role of general relativity in the supernova problem (which appears to be on 
par with other important factors such as the dimensionality and the equation of state) 
both for our understanding of the explosion dynamics as well as for predictions of ob- 
servable signatures. 



1. Introduction 

A core-collapse supernova occurs at the end of the life of a massive star when it has 
exhausted its nuclear fuel at its centre. The nuclear ashes form an iron core which 
undergoes a gravitational collapse once it reaches a mass of roughly \AMq. When 
the core reaches supranuclear densities, the collapse is halted due the stiffening of the 
equation of state, and a shock wave is launched as the core rebounds. The initial kinetic 
energy of the shock is quickly spent for disintegrating the nuclei in the infalling material 
material into free nucleons, and the shock is further weakened by the rapid emission 
of several 10^^ erg in neutrinos from the post-shock region during the neutronization 
burst. The stalled shock typically hovers at a radius of ~ 100 . . . 200 km - possibly for 
hundreds of milliseconds - until it is revived, and expels the envelope of the progenitor. 

The nature of the explosion mechanism has been a subject of intense research in 
computational astrophysics for several decades. Several mechanisms have been pro- 
posed to explain the revival of the shock, the most prominent one being the "delayed 
neutrino-driven mechanism" dBethe & Wilsoniil985. : ,Wilsoa,1985.) . which relies on the 
deposition of energy by neutrinos in the "gain layer" behind the shock to re-energize the 
shock. Except for a special class of progenitors with an 0-Ne-Mg core (iKitaura et all 
l2006h . this mechanism only works in concert with multi-dimensional hydrodynamic ef- 
fects such as convec t ion and the so-called " standing accretion sh ock instability" (SASI, 
Blondin et all l2003l : iFoglizzo et"aLl l2006l : lOhnishi et all l2006l) . which both increase 



the efficiency of neutrino heating in the gain reg ion. W hile some of the most ambi- 
tious supernova simulations in axisym metry (Bura s et"al1 l2b06a: Ma rek & Janka 120091 : 



Bruenn et al. I I2OIOI : Suwa et al.ll2010h indicate that the neutrino-driven mechanism may 



indeed work, there are also alternative scenarios such a magnetohydrodynamically- 
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driven explosions dBisnovatvi-Kogan 19701: Akivama et alj|2003l : iBurrows et al.ll2007l). 
acoustically-powered supernovae ("Burrows et al.ii2006 : see however'Weinberg & Ouat aerj 
2008) for criticism), and explosions triggered by a QCD phase transition (.Sagert et all 



2009D. 



From a computational point of view, core-collapse supernovae present a number 
of challenges: Not only are they inherently multi-dimensional phenomena due to the 
operation of hydrodynamical instabilities, but in order to accurately capture the crucial 
effects of neutrino cooling and heating in the optically thick and thin regimes, the prob- 
lem needs to be treated within the framework of kinetic theory. As the direct solution of 
the full Boltzmann equation for neutrinos is currently only feasible under the assump- 
tion of spherical symmetry dYamada et aLlll999l : iLiebendorfer et alj|2004 ). a variety of 
different approximation strategies for the neutrino transport are used in the most so- 
phisticated multi-dimensional supernova simulation s, including "ray-by -ray" transport 
combined with variable Eddingt on factor techniques (IBuras et al. 2006b) or flux-limited 
multi-group diffusion schemes dBruenn et al.ll2010l) . multi-angle Boltzmann transport 
without energy- bin cou pling (OttetaL I l2008h . two-moment schemes with an analytic 



clos ure (lOb ergaulinger & J ankal 1201 ll) . and the isotropic diffusion source approxima- 
(iLieber 



tion (ILieben dorfer et al.„2009ll 

Moreover, general relativity (GR) plays a major role in the supernova problem 
due to the compactness of the proto-neutron star (GM/Rc^ ~ 0.1 .. . 0.2) and the oc- 
currence of high velocities (~ 0.3c). A general relativistic treatment is also required 
for precise predictions of the gravitational wave signal from core-collapse supernovae. 
However, relativistic supernova simulations with u p-to-date neutrino t ransport have 



However, relativistic supernova simulations witn u p-to-aate neutrino ttansport nave 
long been limited to the case of spherical symmetry ("Brue nn et al.ll200il : lYamada et al. 
[1999; Liebendorfer et al. 2004). In order to study the effects of GR in the context of 
multi-dimensional models, and in particular its i mpact on the gravit ational wave emis- 
sion from the supernovae core, we have recently (^Miiller et al. 20101) introduced a gen- 
eralization of the ray-by - ray-plus variable Eddington factor method (Rampp & Jan^ 
2OO2I : iBuras et aLll2006br) . In this paper, we present results from the first generation 



of axisymmetric relativistic supernova simulations using our Vertex-CoCoNuT code 
and summarize the results of a comparison with the Newtonian case as well as with the 
pseudo-Newtonian "effective potential" approximation. 



2. Relativistic Variable Eddington Factor Method 



In our approach to neutrino transport in core-collap se supernovae , we so lve the equa- 
tions of GR hydrodyna rnics in the formulation of iBanvuls et all (119970 and use the 
xCFC approximation of ICordero-Carrion et al. for the space-time metric, which 

is particularly suitable for our purpose b ecause of its excellent stabilit y properties and 
high accuracy in the core-collapse case (ICordero-Carrion et al. 2009). The hydr ody- 
namics solver is based on the CoCoNuT code (^Dimmelmeier et al. 2002. '20051). an 



implementation of a HRSC scheme with PPM reconstruction ( Colella & Woodward! 
[l984) an d second-o rder Runge-Kutta time-stepping. Different from the original code 
(IMuller et al.ll201(il) . we rely on an improved scheme for mai ntaining tota l energy con- 
servation and employ the relativistic HLLC Riemann solver oflMig none & Bodo. (20051) 
to resolve contact discontinuities in the convective post-shock region more accurately. 

In order to capture the effects of neutrino heating and cooling the equations of 
hydrodynamics need to be coupled with a kinetic equation for the neutrino distribution 
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Figure 1. Left: Evolution of the average shock radius for simulations of a 15Mq 
progenitor with general relativistic hydrodynamics and neutrino transport (black), 
with an effective potential (red), in the purely Newtonian approximation (blue), and 
with a simplified set of neutrino interaction rates in the GR case (light brown). Right: 
Electron antineutrino luminosities (solid) and mean energies (dashed) for the four 
different cases. 



functi on f. Following the approach already used for the Newtonian version the Vertex 
code (|Rampp & Jankall2002l : iBuras et alJl2006bh . we simplify the full six-dimensional 
phase-space problem by considering the first two angular moments J, H,K,L, . . . of 
the energy-dependent radiation intensity and by requiring / to be axially symmetric in 
momentum space (but not in real space) around the unit radius vector (the ray-by-ray 
approximation). The transport problem can thus be reduced to conservation equations 
with source terms for neutrino number, energy and momentum; e.g., the equation for 
the zeroth moment J reads 



d^lyW{J + VrH) d_ 
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in the adopted gauge, with a denoting the lapse function, the conformal factor, yS,- the 
radial shift vector, y the determinant of the three-metric, Vr the radial velocity, W the 
Lorentz factor, and C^'^^ the zeroth moment of the collision integral. A similar equation 
is solved for H, and the higher moments K and L that are required to close the system 
are provided by means of a formal solution of a simplified "model" Boltzmann equa- 
tion. It should be noted that we fully retain the energy-dependence of the moment equa- 
tions, and fully include Doppler shift, gravitational redshift an energy redistribution by 
inelastic scattering of neutrinos off nucleons, nuclei, electrons, and other neutrinos. 



3. Results 



Relativistic supernova si mulations have been conducted fo r two different non-rotatin g 
progenitors with 11.2Mq (IWooslev et al.ll2002D and ISMq dWooslev & Weavej[l995l) . 
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For the ISM© progenitor, additional models were computed to allow for a comparison 
with the purely Newtonian case and the "effective potential" approach /Mare k et~aL 



which mimics certain strong-field effects by modifying the Newtonian gravita- 



tional potential. Moreover, we also considered the case of slightly simplified neutrino 
interaction rates, neglecting the effect of recoil, high-density correlations and weak 
magnetism in neutrino-nucleon reactions and ignoring reactions between different neu- 
trino flavours. As the ISM© progenitor has already prov ed to be a marginal c ase in 
earlier studies relying on the effective potential approach dMarek & Jankal 12009), It IS 
ideally suited to illustrate the impact of slightly different heating conditions depending 
on the treatment of GR. However, the relativistic 

Among the four ISM© simulations, only the GR model with the best currently 
available set of neutrino opacities for our code develops an explosion around 400 ms 
after bounce (left panel of Fig. [T|l. The more optimistic evolution of the GR model com- 
pared to the Newtonian run has been traced to slightly higher surface temperatures of 
the more compact proto-neutron star, which result in higher neutrino luminosities and 
mean energies (right panel of Fig. [T]) and hence allow for more effective heating. To a 
lesser extent, this effect is also present in the effective potential run, but here the faster 
advection (i.e. shorter exposure time) of the accreted material through the smaller gain 
region around the compact neutron star compensates for the increase in the local heat- 
ing rate. On the other hand, the enhancement of the neutrino heating in the GR model 
is strong enough to overcome this competing effect and to shift the balance between 
neutrino heating and advection in the gain region far enough to achieve favourable con- 
ditions for the development of an explosion. We emphasize, however, that despite such 
differences, the effective potential approximation provides a remarkable improvement 
over the purely Newtonian treatment also in multi-dimensional supernova models. 

Incidentally, we also find that the heating c onditions depend qu ite sensitively on 



the neutrino microphysics (as already noted by 'Rampp e t all |2002| and iBruenn et al. 



I2OIQ) as the GR run with simplified interactions rates fails to develop an explosion. The 
more optimistic evolution of the model with the improved rates stems primarily from 
the reduction of the scattering cross-section on nucleons due to weak magnetism 
and nucleon correlations (Horow itz 2002). which helps to enhance luminosities and 
mean energies (right panel of Fig. [1). 

The treatment of GR also turns out to be a crucial factor for the prediction of 
the gravitational wave signal. While the wave signal for our relativistic explosion 
models qualitatively shows the typical features known from (pseudo-)Newtonian stud- 
ies w ith sophisticate d or si mplified neutrino transport dKotake et al.ll2007l : iMarek et"aL 



20091 : iMurphy et al1l2009l : IVakunin et al. .2010) with distinct phases of gravitational 



wave emission from prompt post-shock convection, hot-bubble convection, enhanced 
SASI sloshing motions, and asymmetric shock expansion (left panel of Fig. O, the 
signal spectrum is rather sensitive to GR effects (right panel of Fig. |2l). In the purely 
Newtonian case, the integrated signal (which is dominated by the contribution from 
hot-bubble convection) from the ISM© progenitor peaks at distinctly lower frequencies 
around ~ 500 Hz than in the GR case (~ 900 Hz), which is a consequence of a lower 
Brunt- Vaisala frequency and hence a less rapid deceleration of convective plumes in 
the stably stratified cooling region above the surface of a more extended proto -neutron 
star surface with lower surface gravity (cp. the analysis of i Murphy et al.ll2009 ). On the 
other hand, the median frequency is somewhat overestimated by the effective potential 
approach (~ 1 100 Hz) because the lower neutrino luminosities and mean energies result 
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Figure 2. Left: Matter (black) and neutrino (light brown) gravitational wave sig- 
nals for the general relativistic 15Mq explosion model. Right: Gravitational wave 
spectrum for the first 0.5 s of the post-bounce evolution of a ISM© progenitor for 
simulations depending on the treatment of gravity (black: GR hydro, red: Newto- 
nian hydro + effective potential, blue: purely Newtonian). 



in a steeper density stratification in the cooling region, which in turn leads to a more 
abrupt braking of convective bubbles. The gravitational wave signal is thus at least as 
sensitive to GR effects as to other physical key parameters such as the equation of state 
dMarek et al.ll2009h . 



4. Conclusions and Outlook 



Multi-dimensional general relativistic simulations of core-collapse with a sophisticated 
treatment of the microphysics and the neutrino transport on par with the best currently 
available Newtonian models have only recently become possible, but the first available 
results presented here in this paper already serve to underscore the importance of gen- 
eral relativity in the supernova problem. For a ISM© progenitor, we found that GR 
somewhat improves the heating conditions compared to models computed in the New- 
tonian and the effective potential approximations, which, unlike the GR model, fail to 
explode. The gravitational wave spectra are also considerably changed by GR effects, 
which shift the typical frequency of the time-integrated signal upward by ~ 80% com- 
pared to the purely Newtonian case. We therefore conclude that an accurate treatment of 
GR effects may be no less relevant for a better understanding of the neutrino-driven ex- 
plosion mechanism and quantitative predictions of the signals from core-collapse super- 
novae than other (undoubtedly important) key factors that have recently been discussed 



novae than otner (undoubtedty important) key factors tnat nave recentty been aiscussea 
such as dimensionahty issues dNordhaus et al.ll2O10l : iHanke et al]|201ll : iTakiwaki et al- 



io 111) and the nuclear equation of state. 



Acknowledgments. This work has been supported by the Deutsche Forschungsge- 
meinschaft through the Transregional Collaborative Research Centres SFB/TR27 "Neu- 
trinos and Beyond" and SFB/TR7 "Gravitational Wave Astronomy" and the Cluster of 
Excellence EXC-153 "Origin and Structure of the Universe" (http://www.universe-cluster. 
The simulations discussed here have been performed on the IBM Power6 p575 system 
at the Rechenzentrum Garching (RZG) and the NEC SX-8 at the HLRS Stuttgart. 



6 



MuUer et al. 



References 



Akiyama, S., Wheeler, J. C, Meier, D. L., & Lichtenstadt, I. 2003, ApJ, 584, 954 
Banyuls, E, Font, J. A., Ibanez, J. M. A., Marti, J. M. A., & Miralles, J. A. 1997, ApJ, 476, 221 
Bethe, H. A., & Wilson, J. R. 1985, ApJ, 295, 14 
Bisnovatyi-Kogan, G. S. 1970, AZh, 47, 813 

Blondin, J. M., Mezzacappa, A., & DeMaiino, C. 2003, ApJ, 584, 971 

Bruenn, S. W., De Nisco, K. R., & Mezzacappa, A. 2001, ApJ, 560, 326 

Bruenn, S. W., Mezzacappa, A., Hix, W. R., Blondin, J. M., Marronetti, R, Messer, O. E. B., 

Dirk, C. J., & Yoshida, S. 2010, ArXiv e-prints. 1002 . 4914 
Buras, R., Janka, H.-T., Rampp, M., & Kifonidis, K. 2006a, A&A, 457, 281 
Buras, R., Rampp, M., Janka, H.-T., & Kifonidis, K. 2006b, A&A, 447, 1049 
Burrows, A., Dessart, L., Livne, E., Ott, C. D., & Murphy, J. 2007, ApJ, 664, 416 
Burrows, A., Livne, E., Dessart, L., Ott, C. D., & Murphy, J. 2006, ApJ, 640, 878 
Colella, R, & Woodward, R R. 1984, J. Comp. Phys., 54, 174 

Cordero-Carrion, I., Cerda-Duran, P., Dimmelmeier, H., Jaramillo, J. L., Novak, J., & Gour- 

goulhon, E. 2009, Phys.Rev.D, 79, 024017 
Dimmelmeier, H., Font, J. A., & Muller, E. 2002, A&A, 388, 917 

Dimmelmeier, H., Novak, J., Font, J. A., Ibanez, J. M., & Muller, E. 2005, Phys.Rev.D, 71, 
064023:1 

Foghzzo, T., Scheck, L., & Janka, H.-T. 2006, ApJ, 652, 1436 

Hanke, F, Marek, A., Mueller, B., & Janka, H.-T. 201 1, ArXiv e-prints. fTTS8 . 43551 

Horowitz, C. J. 2002, Phys.Rev.D, 65, 043001 

Kitaura, E S., Janka, H.-T., & Hillebrandt, W. 2006, A&A, 450, 345 

Kotake, K., Ohnishi, N., & Yamada, S. 2007, ApJ, 655, 406 

Liebendorfer, M., Messer, O. E. B., Mezzacappa, A., Bruenn, S. W., Cai-dall, C. Y, & Thiele- 
mann, F-K. 2004, ApJS, 150, 263 
Liebendorfer, M., Whitehouse, S. C, & Fischer, T. 2009, ApJ, 698, 1 174 
Marek, A., Dimmelmeier, H., Janka, H.-T, Muller, E., & Buras, R. 2006, A&A, 445, 273 

Marek, A., & Janka, H. 2009, ApJ, 694, 664. 0708 .3372 

Marek, A., Janka, H., & Muller, E. 2009, A&A, 496, 475. 08 08.41361 

Mignone, A., & Bodo, G. 2005, MNRAS, 364, 126 

Muller, B., Janka, H.-T, & Dimmelmeier, H. 2010, ApJS, 189, 104. 11001 . 48"4T1 

Murphy, J. W., Ott, C. D., & Burrows, A. 2009, ApJ, 707, 1 173. 090 7.47621 

Nordhaus, J., Buitows, A., Almgren, A., & Bell, J. 2010, ApJ, 720, 694. 11006.3792] 

ObergauHnger, M., & Janka, H.-T. 2011, ArXiv e-prints. 1101 . 1198 

Ohnishi, N., Kotake, K., & Yamada, S. 2006, ApJ, 641, 1018 

Ott, C. D., Buitows, A., Dessait, L., & Livne, E. 2008, ApJ, 685, 1069 

Rampp, M., Buras, R., Janka, H., & Raffelt, G. 2002, in Nuclear Astrophysics, edited by 

W. Hillebrandt & E. Muller, 1 19. arXiv : a stro-ph/0203493| 
Rampp, M., & Janka, H.-T. 2002, A&A, 396, 361 

Sagert, I., Fischer, T, Hempel, M., Pagliara, G., SchafFner-Bielich, J., Mezzacappa, A., Thiele- 
mann, F-K., & Liebendorfer, M. 2009, Physical Review Letters, 102, 08 1 101 : 1 

Suwa, Y, Kotake, K., Takiwa ki, T., Whitehouse, S. C, Liebendorfer, M., & Sato, K. 2010, 
PASJ, 62,L49h-. 0912.TT57 

Takiwaki, T., Kotake, K., & Suwa, Y. 201 1, ArXiv e-prints.[ IT08 .3989 1 

Weinberg, N. N., & Quataert, E. 2008, MNRAS, 387, L64. 0802 . 1522 

Wilson, J. R. 1985, in Numerical Astrophysics, edited by J. M. Centrella, J. M. Leblanc, & R. L. 
Bowers, 422 

Woosley, S. E., Heger, A., & Weaver, T. A. 2002, Rev. Mod. Phys., 74, 1015 
Woosley, S. E., & Weaver, T A. 1995, ApJS, 101, 181 

Yakunin, K. N., Marronetti, P., Mezzacappa, A., Bruenn, S. W., Lee, C, Chertkow, M. A., Hix, 
W. R., Blondin, J. M., Lentz, E. J., BronsonMesser, O. E., & Yoshida, S. 2010, Classical 
and Quantum Gravity, 27, 194005. 10 05.0779 

Yamada, S., Janka, H.-T., & Suzuki, H. 1999, A&A, 344, 533 



